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Massless Dirac fermions occur as low-energy modes in several quasi-two-dimensional condensed 
matter systems such as graphene, the surface of bulk topological insulators, and in layered organic 
semiconductors. When the rotational symmetry in such systems is reduced either by an in-plane 
electric field or an intrinsic tilt of the Dirac cones, the allowed dipolar optical transitions evolve 
from a few selected transitions into a wide fan of interband transitions. We show that the Lorentz 
covariance of the low-energy carriers allows for a concise analysis of the emerging magneto-optical 
properties. We predict that infrared absorption spectra yield quantitative information on the tilted 
Dirac cone structure in organic compounds such as a-(BEDT-TTF) 2 l 3 . 

PACS numbers: 73.22.Pr,73.61.Ph,78.20.Ls 


I. INTRODUCTION 


The low-energy carriers in graphene imitate the rel¬ 
ativistic physics of massless Dirac fermions with the 
Fermi velocity v « c/300 replacing the speed of light 
Cr^ Their peculiar Landau level (LL) structure provided 
the first evidence for the successful fabrication of mono- 
layer graphene.— Experimental progress in the last decade 
has allowed, for example, the direct measurement of the 
Dirac cones in photoemission experiments,— the observa¬ 
tion of LLs in magneto-optical^ and scanning-tunneling 
experiments^ and the detection of Klein tunneling.— Fur¬ 
ther predictions follow from the Lorentz covariance of 
low-energy electrons in graphene, notably the collapse 
of Landau levels subjected to in-plane electric and per¬ 
pendicular magnetic fields.— Here, conventional Landau 
states exist only as long as a Lorentz boost is possible to a 
reference frame where the external electric field vanishes; 
both the spectrum and the eigenstates follow by Lorentz 
covariance. Similar ideas apply to massless Dirac surfaces 
states of bulk topological insulators.— On the theory side, 
exploiting the Lorentz covariance of the massless Dirac 
equation often yields very concise derivations (compare, 
e.g.. Refs. 0 and [H, with the additional advantage of a 
unified understanding of a class of high-energy and low- 
energy phanomena. 

More recently, the quasi-2D organic compound a- 
(BEDT-TTF) 2 l 3 has been shown to host Dirac cones at 
low energie o^^d^ under pressurei^ or uniaxial strain—^ 
Other quasi-2D organic materials that may share the 
same property, e.g., 0 — (BEDT-TTF) 2 l 3 ^‘^ and a- 
(BEDT-TTF) 2 NH 4 Hg(SCN) 4 ,— still need further inves¬ 
tigation. In contrast to graphene, the Dirac cones in 
these materials are tilted (see Fig. [T]). We will refer to 
such compounds as 21) Weyl materials. In a magnetic 
field, this tilt plays the role of an effective in-plane electric 




FIG. 1. (Color online) Organic compound a-(BFDT- 
TTF) 2 l 3 . (a) Molecule (BEDT-TTF) that consists of car¬ 
bon (C) and sulfur (S) atoms, (b) Arrangement of the four 
molecules (A, A’, B, and C) in the unit cell of a single layer 
(a6-plane). (c) Sketch of the low-energy band structure in 
the vicinity of the Fermi level that has the form of two tilted 
Dirac cones. We set up the coordinate system so that the 
cones are tilted in the direction. 


fieldii and can be treated, as we show here, in a covariant 
manner. Moreover, the covariance of Dirac electrons in 
graphene and in a-(BEDT-TTF) 2 I 3 yields another rel¬ 
ativistic effect that consists of highly unusual magneto¬ 
optical transitions. Whereas in the frame of reference, 
where the electric field (or the tilt in 2D Weyl materi¬ 
als) vanishes, one obtains the usual dipolar transitions 
between the LLs ±n and ±(n ± 1),— we predict that a 
Lorentz boost back to the laboratory frame gives rise to 
an amazing multiplication of these optical lines into a fan 
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of experimentally observable interband transitions. Fi¬ 
nally, we show that magneto-optical measurements may 
yield information on the tilt of the Dirac cones in 2D 
Weyl materials. 

This article is structured as follows. In Sec. El we re¬ 
view the low-energy band structure of the materials we 
study together with the relevant single-particle (Landau) 
states. In Sec. lIIIl we study the matrix elements of the in¬ 
teraction of the electrons with electromagnetic radiation. 
The optical absorption per layer is analyzed in Sec. IIVI 
where we also discuss how information can be gathered 
on the tilted Dirac cone structure of 2D Weyl materials 
from magneto-optics. In Sec.|V]we conclude and discuss 
the experimental connections. Details of the calculation 
are delegated to the Appendices. 


II. MATERIALS AND MODELS 

Let us first discuss graphene electrons in a perpendic¬ 
ular magnetic field B = and in-plane electric held 

E = —Eyey. The low-energy carriers are described by 
the Hamiltoniaiii 

Hg = v{p + eA) ■ a - eE^yl, (1) 

where v is the Fermi velocity of graphene, the Pauli ma¬ 
trices are denoted by at with i G {x,y,z} and cr = 
(cTj;, ( 7 j,). Hereafter, we use the Landau gauge A = 
—yB±ex and e > 0. With the choice of = {vt, x, y) in 
the laboratory frame 'R- one may introduce the 2D mass¬ 
less Dirac equation 


Ihv^^ j = 0, (2) 

where the representation 70 = Uz, 71 = —(Jz<Jx, 72 = 
—(Jz(Ty has been used, and = {E^^y/v, —yBj_, 0). The 
electric held can be eliminated by a Lorentz boost in the 
^-direction from the laboratory frame TZ to {vf, x', y') = 
A(vt,x,y) = {j{vDt+l3x),j{vDt+x)^y) to another refer¬ 
ence frame TZ', as long as the drift velocity vd = E^\/B± 
is smaller than the Fermi velocity, vd < v. The Lorentz 
boost is thus parametrized by the rapidity tanh 9 = (3 = 
Vd/v and the Lorentz factor 7 = l/\/l — The Hamil¬ 
tonian in TZ' then reads 


H'g = v(p' -f eA') • 0 -, (3) 

and covariance requires the two-spinor to transform as 
'^'{vt',x',y') = S{A)'^{vt,x,y), 
with S'(A) = exp(|cr,j). 

Notice that the simple form of Hq allows for 
a straightforward diagonalization in TZ' by standard 
techniquesThe sought eigenstates of Hq are then 


obtained by a transformation back to TZ^ 


'^\,n-k{x,y) 


+A 


sinh I 
— cosh I 

^ cosh I 
^ — sinh I 


(p\,n-k{x,y) 


^X,n—l;k {x^ y) 


( 4 ) 


for n > I and 

4'n=o;fe(a;,y) = ^ K=o-,k{x,y) (5) 

for n = 0, where A = ± denotes LLs of positive and neg¬ 
ative energy, respectively, and k is the wavevector in the 
x-direction related to the conserved momentum compo¬ 
nent px = hk. Here we have used the functions 

(l)\,n-k{x,y) = ^-J= (7?A,n) (6) 

V27r ^y/TT^i2'^n\ 

V\,n = (y - kf) - XlSV^, (7) 

where i = yJh/eB is the magnetic length, Hn is a Her- 
mite polynomial, and i5-normalization has been used for 
the wave functions. The covariant character is apparent 
here in the argument of the wavefunctions, where 
one notices the contraction factor Xjyfp. The second 
term reflects the admixing of the time-like component, 
A-\/2n, proportional to the LLs, 

Q XHv f —— . , 

^A.n;fc = 2n - hvDk, (8) 

which exhibit a collapse in the ud —^ u (7 —^ 00 ) limitji'i^ 
The last term in Eq. ([S]) stands for the dependence of 
the single-particle energies on the location of the guiding 
center if an in-plane electric field is present. This term is 
also present in the case of nonrelativistic electrons, but 
the LL spacing of the latter is insensitive to the inplane 
electric field, in contrast to the decreased LL spacing of 
relativistic electrons apparent in the factor Iin Eq. 
( 0 ). 

The low-energy carriers in 2D Weyl materials are de¬ 
scribed by the minimal Weyl Hamiltonianfiii^SiiSi which 
allows for different velocities Vx, Vy in two directions, and 
a shift of the centers of the equienergy ellipses in a generic 
tilt direction. The anisotropy Vx yf Vy can be eliminated 
by rescaling y' = (vx/vy)y, and rotating the direction of 
the tilt into the a;-direction. The details of this procedure 
are presented in Appendix We obtain 

Hw = Vx{p + eA) ■ a-+ vq{px E eAx)l. ( 9 ) 

Here vq is a velocity parameter that is related to the band 
structure, c.f. Eq. (RtI) . 

In the Landau gauge the noncovariant term VqPx^ in 
Eq. (|n]) is harmless as Px only extracts a quantum num¬ 
ber. The Hamiltonian ([H]) can be rewritten in the form 
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i?w = + foPxl, where the last term is diagonal 

both with respect to the sublattice cr and in the wave 
vector k. The energy levels can therefore be written as 


e 


w 

X,n;k 


cov. 

^A,n;/c 


+ hvok, 


where are the solutions of the Hamiltonian’s co¬ 
variant part which has the same form as that 

in Eq. O- In the Landau gauge, the term ev^A^l = 
—ev^Bx^yl. can be viewed as generated by an effective 
pseudo-electric field Ey® = and Vq assumes the 

role of the drift velocity. The wave functions are ob¬ 
tained by methods that are completely analogous to the 
graphene case, using tanh0 = (5 = vq/v m. the Lorentz 
boost. Then the LL spectrum becomes^SiSi 


w 

^X.n-.k 


Xhv 

73 / 2 ^ 




( 10 ) 


where an average velocity v = JVxVy has been intro¬ 
duced. Comparing with Eq. the guiding-center- 
dependent term is absent, because there is no inplane 
electric field. 

Both in graphene and in 2D Weyl materials, Dirac 
fermions occur in two cones related by time-reversal sym¬ 
metry. Whereas electrons in both valleys are subjected 
to the same external electric field (drift velocity), the tilt 
in Weyl materials is opposite in the two valleys and so is 
the implied effective electric field. 



FIG. 2. (Color online) (a,b) The matrix elements P 

for graphene in crossed B± and i5|| fields, for right circularly 

polarized light. (c,d) for a 2D Weyl material, for el¬ 

liptical polarization that appears as right circularly polarized 
in the rotated and rescaled system. (a,c) /3 = 0.2 and (b,d) 
P = 0.8. Inset of panel (a): the same as in (a) for = 0 
limit. In panel (b) the dotted line indicates the matrix el¬ 
ements corresponding to n' = n > 1,A' = —A transitions, 
see Eq. m-, blue (red) curves indicate equi-energy lines with 
minimal (maximal) value of the matrix elements. 


III. INTERACTION WITH 
ELECTRO-MAGNETIC RADIATION 


Let the 2D sample be located in the z = 0 plane. We 
consider the interaction with a classical radiation field 


A‘'*^®(r,t) =Re 


Eq 

UJ 


{a&x + TBy) 


( 11 ) 


where \a\^ + \t\^ = I and uj is fixed; A''’^®(r, t) is added to 
the vector potential in Eqs. ([T]) and dH). Notice that the 
different single-particle Hamiltonians for graphene and 
2D Weyl materials yield slightly different perturbation 
operators 


^5H = ev—{aux + TUy) (12) 

and 

5H = evx— f atJx + roy H—^51 ) , (13) 

w V Vx J 

respectively. The parameters a and r determine the light 
polarization, e.g., a = l/\/2 and r = for circular 

polarization O / O. In the case of 2D Weyl materials, a 
and T are also affected by the rescaling of the coordinates 
(c.f. Appendix El) : the use of tildes for a and f refer to 
this change. 


For the electron-light interaction, we will need the ma¬ 
trix elements of 

^ - k') (“/“'cf;"'), (14) 

where the dimensionless factors ^ and n of 
these matrix elements are 


x,u =- 7 ; K,u' al3 + XXF^L{ j_^al3 


1 

2 1 


I \ / TTiA.n.A', n'/ I • —1\ I \ TTiA.n.A', n' / • —1\ 

+ [a + lT-f )+AE„12 ; {a - ITJ ) 


y, 


(15) 


and 


w. 




Qx,n =7^ (7 


X,n^X' ,n' / —1 




a — if) 
(7~^ 


a - 


-I-if) . (16) 


The functions •p are defined in terms of associated 

Laguerre polynomials L™, as 


rpX,s,X' ,p 



^xf = pppp - A\/s), (17) 
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for q> r, otherwise 'P = 

We show P in Fig. [5] for graphene and 2D 

Weyl materials for /3 = 0.2 and 0.8, assuming right cir¬ 
cular polarization and a special choice of elliptical polar¬ 
ization, respectively. The latter choice is determined by 
the asymmetric velocities and corresponds to a circular 
polarization after rescaling. The opposite polarization is 
obtained via the relation 

and both \{Q'^^)xn P I(*3 ^)a n" 1^ independent 
of the sign of f3. The results for both systems are al¬ 
most identical, apart from small differences due to the 
slightly different perturbation operators in Eq. (HI and 
m- Most saliently, we obtain novel transitions beyond 
the usual dipolar ones, oc dn',n±i that arise for vanishing 
i?ll or no tilt, respectively, and which apply therefore only 
in the reference frame TV, where rotation symmetry is re¬ 
stored. Alternatively, one may thus view the occurrence 
of additional transitions as due to the broken rotation 
symmetry caused by the Lorentz boost back to the lab¬ 
oratory frame TZ. The weight of intraband transitions 
that change the LL index n by more than one unit grows 
monotonically with P whereas the interband transitions 
show a complex fan structure with an arc pattern. The 
fan, whose opening angle increases with /3, is bounded 
by large values of the perturbation matrix element that 
are not affected by the arc pattern [c.f. also the inset of 
Fig.Hb)]. 



FIG. 3. (Color online) The matrix element as a 

continuous function of n (green) for /3 = 0.8. The black dots 
indicate integer values of n. Inset: the same along the lines 
\' = —A and n' = n -\- c, where c = 1, 2, 3. 

For further insight into the origin of the observed arc 
pattern, consider the case n'=n>I,A' = —A for large 
p. One finds that 

( 19 ) 



FIG. 4. (Golor online) (a) Absorption coefficient of graphene 
in perpendicular magnetic and in-plane electric fields for right 
circularly polarized light. Curves are shifted for clarity, (b) 
Comparison of iabs for right and left circularly polarized light 
for P = 0.7. Inset: The total labs, bnner due to the inner 
part of the interband fan, and iouter due to the outer part, for 
graphene /3 = 0.7. 


where the associated Laguerre polynomial L™ can be 
generalized to continuous n in terms of the confluent 
hypergeometric function. At fixed /3, this function has 
decaying oscillations in n. For P = 0.8, this function is 
plotted in Fig. [3Ka), where we have used n as a continu¬ 
ous variable. The apparent period 5n is more than, but 
close to, unity, and this incongruity with integer values 
yields a beating with a distance An ~ 10 between consec¬ 
utive zeros that gives rise to the observed arc pattern. By 
Nyquist’s sampling theorem, this effect survives as long 
as 6n < 2 (which corresponds to /3 > 0.4). Nearby transi¬ 
tions such as (A',n') = (—A,n-|-1), (A',n') = (—A,n-|-2), 
etc., can be analyzed similarly. Remarkably, their beat¬ 
ings are of almost identical period and are in phase (Fig. [3] 
inset). While the latter phase condition does not hold for 
small P in general, we will see below that two characteris¬ 
tic periods can be extracted, one of which is continuously 
connected to the period at large p. 
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IV. MAGNETO-OPTICAL ABSORPTION 

By Fermi’s Golden Rule;^ the transition rates per unit 
area A are 

pcmi/abs _ 1 ^TT Oc ttX’ ^k' ^ 

,n'-k' ~ ^^\,n,k 

X S{ey,n' - e\.n ± fia;)nF(eA,n)[l - npiey^n')], (20) 

where np^e) is the Fermi function. Integrating over all 
states, the total transition rate per unit area is 

R= Y. fdk fdk' - Rl^l';k'^X,n;k] 

( 21 ) 

This can be normalized using the absorbed/emitted en¬ 
ergies and the magnitude of the Poynting vector S = 
E^/2fj,oC, 


•G/w 

^abs 


(w) 


Rhuj 


( 22 ) 


which is dimensionless. (The celebrated 2.3% per sheet 
absorbtion of pristine graphene^ employs this defini¬ 
tion.) We introduce a Lorentzian broadening F for ne¬ 
glected effects (temperature, impurities, phonons, etc). 
Taking spin and valley into account, one obtains 


G/w 

^abs 


(^) 


2v^hfj,oce^ 

P 


E 

■n . \' 


nF{e\,n) - 


X 


G/WQX',n' 


21 r 

TT r2 -I- (eA',„' - e\.n 


nF{e\',n') 

~ €A,n 


hwY 


(23) 


Inspecting the graphene data by Sadowski et al^ we will 
use F = 0Ahv/£ as a rough estimation. Eq. (TTSl) im¬ 
plies that the absorption for opposite polarizations are 
connected 


■G/W.O 

^abs 


(w,Ai) 


■G/W.O/ N 

*abs (w,-M), 


(24) 


where /i is the chemical potential. 

The absorption of graphene in perpendicular megnetic 
and in-plane electric field is plotted as a function of fre¬ 
quency in Fig. UKa) for several values of /3. Whereas 
for low frequencies one finds well-pronounced lines corre¬ 
sponding to individual transitions, in the large frequency 
domain the lines of many allowed transitions coalesce. 
The resulting supermodulation stems from the arc pat¬ 
tern in Fig. [21 where the arcs connect different transi¬ 
tions of nearby energies, see the lines connecting min¬ 
ima/maxima in Fig.|2](b) for illustration. In Fig.ld^b) we 
compare the absorption for the two circular light polar¬ 
izations. The low-frequency absorption, where only few 
LL transitions contribute, depends strongly on polariza¬ 
tion. In contrast, the absorption becomes polarization- 
independent at high frequencies because of the coales¬ 
cence of a large number of transitions. However, the ef¬ 
fect of the above-mentioned supermodulation is retained. 


Finally, in the inset of Fig.|4Kb) we separate the contribu¬ 
tions of the patterned inner wedge of the fan of interband 
transitions and of the bounding lines to *abs(w); we find 
that the latter merely yields a constant background. 



FIG. 5. (Color online) (a) The one-sided Fourier-transformed 
rescaled absorption spectrum J^[iabs(ta')](t) for linearly po¬ 
larized light exhibits two characteristic peaks in the range 
/3 < 0.3. We show the case of 2D Weyl materials; the graphene 
case is similar, (b) For, /3 > 0.4, only a single characteristic 
peak is discernible. Inset: the dependence of the character¬ 
istic peak on the direction of linear polarization at d = 0.6, 
relative to the tilt direction. 

All features in the optical spectra follow the o'- 

^/n scaling of massless 2D Dirac fermions. Therefore, we 
consider the spectra as a function of uj' = and seek 
the effect of the high-frequency regular pattern in their 
one-sided Fourier transform (or Laplace transform with 
an imaginary argument), 

'iabs(a; = v^). (25) 

A characteristic peak emerges, as shown in the inset of 
Fig. IS] for 2D Weyl materials, even for linearly polarized 
light. A particular polarization direction exists for which 
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the visibility is optimal. For graphene, this is the direc¬ 
tion of the electric field; for 2D Weyl materials, it is the 
tilt direction. For small /? < 0.35, a satellite peak ap¬ 
pears [c.f. Fig. [5l[a) for the case of 2D Weyl materials], 
but the main peak remains clearly visible. 



P 


FIG. 6. (Color online) (a) The location of the main peak(s) 
in the one-sided Fourier-transform of the rescaled absorption 
spectrum J^[iabs(w0](^) for linearly polarized light as a func¬ 
tion of the tilt parameter j3 for a 2D Weyl material. Error 
bars indicate the FWHM assuming a Lorentzian broadening, 
(b) The same data plotted in units of where the 

angle-averaged Fermi velocity v* = yJVxVy{l — is di¬ 

rectly related to the density of states that appears in several 
experiments. 

Figure [5] shows the relation between the tilt parame¬ 
ter (3 and the location of the two peaks in the one-sided 
Fourier-transform of the rescaled spectrum in Eq. (1251) . 
Here, we restrict our attention to 2D Weyl materials. Ex¬ 
tracting /3 is straightforward in the moderate-tilt limit, 
/3 < 0.35, via the splitting of the two peaks in J'[zabs(w')]. 
In the large-tilt region /3 > 0.35, however, there is only 
one peak. Its location is informative, but for quantita¬ 
tive analysis we also need to know the Eermi velocity. 
Then it is more useful to plot the location of the peak 
in units of where the effective Eermi velocity 


V* = y/VxVyh'^^'^ = ^ takes into account 

both the anisotropy and the relativistic reduction of the 
energy scale. This velocity v* occurs directly in the den¬ 
sity of states^i^ or in the LL spacing in Eq. (fTUl) . Thus 
we can obtain the tilt parameter only if we utilize in¬ 
dependent information on v* from thermodynamics or 
low-frequency magneto-optics. 

V. CONCLUSION 

In conclusion, the pseudo-relativistic nature of elec¬ 
trons in graphene in crossed electric and magnetic fields, 
and in 2D Weyl materials such as q;-(BEDT-TTF) 2 I 3 
manifests itself in unique magneto-optical properties. A 
large number of transitions beyond the usual dipolar ones 
become possible. Their coalescence at high frequencies is 
the fingerprint of the broken rotation invariance in both 
cases either due to the particular tilt direction or that of 
the electric field. 

In graphene, this effect might be observed easier than 
the earlier predicted Landau level collapseji because the 
magneto-optical effect is present in the low-/3 range where 
the change of the Landau level energies is still small. A 
side-gated geometry must be necessary. Here, however, 
the screening of the external electric field by the edge 
states22i2^ might be a complication)^ We note that an 
upper limit to the achievable (3 values in graphene is set 
by the condition that potential energy change over a mag¬ 
netic length should be less then the energy difference be¬ 
tween adjecent LLs to avoid tunneling along the parallel 
electric field. For the n = 0,1 Landau levels, this means 
eE\\^ < ef - = ^(1 - which sets /3 < 0.75. 

In a-(BEDT-TTF) 2 l 3 and potentially other quasi-2D 
organic materials with massless Dirac conesj^^— on the 
other hand, a quantitative analysis of this fingerprint, in 
combination with other information on Fermi velocities 
if the tilt turns out to be large, helps us determine the 
tilting parameters of the massless Dirac cones. We no¬ 
tice finally that the same magneto-optical features should 
appear in strained graphene, where uniaxial strain yields 
tilted Dirac cones^S that are predicted to show a clear 
signature also in Raman spectroscopyi^ 
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Appendix A: Derivation of the minimal Weyl 
Hamiltonian 



by the parameter 

V = \J{v^/vxY + {vl/vyY- (A2) 

Following Refs. [U and [13, we use a rescaled and ro¬ 
tated coordinate system, defined by the transformation 



in terms of the rotation angle 4>. Rescaling py removes the 
anisotropy, and the rotation brings px in the tilt direction 
if we choose 


cosp 


VQVy 


(A5) 


Then the four-parameter Hamiltonian Hy/ can be rewrit¬ 
ten in a simpler form 


FIG. 7. (Color online) The one-sided Fourier transform of 
for the Weyl Hamiltonian for linearly polarized light 
in the tilt direction [Eq. (lASIl ] and /3 = 0.5 for various Landau 
level cutoff values: 25, 50, 100, 150, which corresponds to B = 
12.8,6.4,4.2,3.2 T, respectively. Inset: the original 
function. 


Hyj = VxP ■ a-+ vopxl, (A6) 

where p is measured from the direction of the tilt, and 

vo=r]Vx. (A7) 

The Peierls substitution in Eq. (IA6I) then yields Eq. (El- 


Massless Dirac carriers of a-(BEDT-TTF) 2 l 3 are de¬ 
scribed by the minimal Weyl Hamiltonianii using four 
parameters Ug, Wq , , t'y as follows 


Hw = VxPxCTx + VyPyO-y + {v^Px + VoPy)'^' (^1) 


The inclination of the Dirac cone is determined by the 
combined effect of the tilt and the anisotropy. We char¬ 
acterize the anisotropy of the Dirac cones by the quotient 
Vxjvy. Even if Vx = Vy, the constant energy slices are not 
concentric because (v§,Vq) ^ (0,0). We quantify this tilt 


Appendix B: Determination of the cutoff 

As the number of Landau levels within the range of 
validity of the Weyl Hamiltonian is limited, we examine 
the robustness of our results against the change of the 
Landau level cutoff. Previously, we estimated the num¬ 
ber of available Landau levels is SOO/Hj..— Fig. [3 shows 
the effect of an arbitrary change of the cutoff. The os¬ 
cillatory behavior remains visible in the *^s(w) curves 
even at low cutoffs. The relevant peaks of J^[fabs(w')] 
also survive, though with reduced height. We can see 
that the oscillatory behavior of is present for as few 
as 50 Landau levels, which is a safe choice for reasonable 
magnetic fields. 
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